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SYSTEM AND METHOD OF GENERATING A FINITE ELEMENT MESH 
5 FOR A THREADED FASTENER AND JOINING STRUCTURE 

ASSEMBLY 

BACKGROUND OF THE INVENTION 

10 

1 . Field of the Invention 

The present invention relates generally to 
f5 15 finite element analysis of a threaded . fastener and, 

more specifically, to a system and method of 
s'j generating a finite element mesh for a threaded 

_~ fastener and joining structure assembly. 

il 20 2 . Description of the Related Art 

f«s Fastener design, and in particular the 

design of a threaded fastener, has advanced to a 
state in which computer-aided design techniques are 
frequently incorporated in the development and 

25 analysis of a threaded fastener and a structure 
joined together by the threaded fastener, such as two 
blocks. Computer-aided design is especially 
beneficial in analyzing a stress in a joining 
structure resulting from the load transfer through 

30 contacting threads of the threaded fastener and 
mating surfaces in the joining structure. For 
example, a fatigue failure, such as a crack 1 in a 
thread root, may occur in a threaded portion of the 
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joining structure. The ability to predict a stress 
and potential fatigue failure in a threaded fastener 
or joining structure under an operating load is 
advantageous in the design of the fastener and the 
5 bolted joint . 

It is known to utilize a strain gauge to 
determine the state of stress on a part in which a 
load has been applied. However, it is extremely 

C3 

:= difficult to position a strain gage at the thread 

*jf 10 interface between the threaded fastener and the 

joining structure and still be able to recreate the 
stress. It is also difficult to reproduce the effect 
:J of other loading factors, including thermal loading. 

\2 It is also known to use a classical 

15 analytical technique to identify stresses and areas 
of high stress concentration within a structure. For 
example, the thread and joining structure behavior is 
approximated and the model is empirically corrected 
to fit experimental observations . However, this 
20 technique only evaluates stresses and fatigue in the 
bolt shank, bolt thread stripping, and bolt fatigue, 
and doesn't consider stress in the block threads or 
contacting surfaces . 

Another known analytical technique is 
25 finite element analysis. In the past, a simple model 
was used to reduce model size and complexity. For 
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example, the fastener and the joining structure are 
modeled in two dimensions, and a computer simulation 
is performed to approximate the stresses of the 
contacting threads. However, a two-dimensional 

5 analysis assumes that the bolt is axi-symmetric , and 
helical threads are not axi-symmetric. Thus, there 
is a need in the art for a system and method of 
generating a finite element mesh of a threaded 
fastener and a joining structure for stress analysis 
10 and evaluation purposes, to more accurately predict 
fatigue behavior. 

SUMMARY OF THE INVENTION 

Accordingly, the present invention is a 
15 system and method of generating a finite element mesh 
for a threaded fastener and joining structure 
assembly. The system includes a computer system for 
generating a mesh model of the threaded fastener and 
joining structure assembly by creating nodes and 
20 elements for each non-threaded portion of the 
threaded fastener and joining structure assembly 
using cylindrical coordinates and creating nodes and 
elements for each threaded portion of the threaded 
fastener and joining structure assembly using helical 
25 coordinates . The system also includes a user 

evaluating the mesh model of the threaded fastener 
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and joining structure assembly using finite element 
analysis and evaluating a result of the finite 
element analysis. The system further includes a user 
predicting a stress of the threaded fastener and 
5 joining structure assembly from the evaluation of the 
result of the finite element analysis. 

The method includes the steps of generating 
a mesh model of the threaded fastener and joining 
structure assembly by creating nodes and elements for 

10 each non-threaded portion of the threaded fastener 
and joining structure assembly using cylindrical 
coordinates and creating nodes and elements for each 
threaded portion of the threaded fastener and joining 
structure assembly using helical coordinates. The 

15 method also includes the steps of evaluating the 
mesh model of the threaded fastener and joining 
structure assembly using finite element analysis and 
evaluating a result of the finite element analysis. 
The method further includes the steps of predicting a 

20 stress' of the threaded fastener and joining structure 
assembly from the evaluation of the result of the 
finite element analysis. 

One advantage of the present invention is 
that an improved system and method of generating a 

25 finite element mesh for a threaded fastener and 
joining structure is provided that produces a good 



quality three dimensional mesh while considerably 
reducing design time and related expenses. Another 
advantage of the present invention is that the system 
and method utilizes parametric automated design in 
light of predetermined engineering and manufacturing 
criteria. Yet another advantage of the present 
invention is that the system and method facilitates 
fatigue analysis of the threaded fastener and joining 
structure early in the design process. Still another 
advantage of the present invention is that the system 
and method supports computer-aided engineering 
analysis (CAE) and rapid prototyping. A further 
advantage of the present invention is that the system 
and method enhances flexibility in design, while 
still meeting vehicle timing considerations . Still a 
further advantage of the present invention is that 
the system and method utilizes knowledge, guidelines 
and lessons learned from design, engineering and 
manufacturing experience to design the threaded 
fastener and joining structure to improve the 
quality, reduce development time and reduce the cost 
of such. 

Other objects, features and advantages of 
the present invention will be readily appreciated, as 
the same becomes better understood, after reading the 
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subsequent description taken in conjunction with the 
accompanying drawings . 



BRIEF DESCRIPTION OF THE DRAWINGS 

5 FIG. 1 is a block diagram of a system which 

may be utilized with a method of generating a finite 
element mesh for a threaded fastener and joining 

f=i structure, according to the present invention. 

Iq 

™ FIG. 2 is a cross-sectional view of a 

y * 

Ltf 

.V: 10 threaded fastener and joining structure, according to 

ff% 

%a the present invention. 

FIGS. 3A through 3R are flowcharts of a 
f3 system and method of generating a finite element mesh 

m for a threaded fastener and joining structure, 

?g 15 according to the present invention. 

FIGS. 4A through 4C are examples of a 
finite element mesh for a threaded fastener and 
joining structure, according to the system and method 
of FIGS. 3A through 3R. 
20 FIGS. 5A through 5E are examples of a 

uniform finite element mesh for a threaded fastener 
according to the system and method of FIGS. 3A 
through 3R. 

FIGS. 6A through 6C are examples of a non- 
25 uniform finite element mesh for a threaded fastener 



according to the system and method of FIGS. 3A 
through 3R. 

FIG. 7 is an example of a transition 
profile mesh for a threaded fastener according to the 
system and method of FIGS. 3A through 3R. 

FIGS . 8A through 8C are examples of a 
transition profile mesh for a threaded block 
according to the system and method of FIGS. 3A 
through 3R. 

DESCRIPTION OF THE PREFERRED EMBODIMENT ( S ) 

The design of a fastener and joining 
structure, and in particular the design of a threaded 
fastener and joining structure for use as an assembly 
on a vehicle, is achieved according to the present 
invention with a generic, parametric driven design 
method. Advantageously, this system, and method 

allows flexibility in design of the threaded fastener 
and joining structure and engineering analysis of the 
design in a fraction of the time required using 
conventional design methods, since the design is 
automatically evaluated against rules in a knowledge 
base. Various computer-based tools are integrated 
into a single user interface to achieve this enormous 
time and expense savings, including solid modeling, 
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parametric design, automated studies and a knowledge- 
based engineering library. 

Referring to the drawings and in particular 
FIG. 1, the tools 10 used by a system and method for 
5 generating a finite element mesh for a threaded 
fastener and joining structure, according to the 
present invention, are illustrated graphically. The 
tools 10 include a knowledge-based engineering 
library 12 stored on an electronic storage device 

10 (not shown) . The knowledge-based engineering library 
12 includes design, engineering, and assembly rules 
for a fastener or joining structure. In this 

example, the knowledge -based engineering library 12 
is a database of sub-libraries containing an 

15 electronic representation of various experts' 
knowledge of information relevant to the design of 
the fastener and joining structure. For example, the 
knowledge -based engineering library 12 includes a 
component parts library such as a fastener library 

20 containing a database of various types of fasteners 
that may be available for use with the joining 
structure. It should be appreciated that the fastener 
may be parametrically modeled. The knowledge -based 
engineering library 12 may also provide interactive 

25 access to other web-based libraries. 
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The tools 10 also include a vehicle 
platform library 14 stored on the electronic storage 
device. The vehicle platform library 14 is an 
electrical representation of a vehicle platform or a 
5 portion thereof. For example, the vehicle platform 
library 14 may include a model of a particular 
vehicle body design, or a portion thereof such as an 
engine mount or engine block. It should be 

appreciated that the vehicle platform library 14 may 

10 be a sub-library within the knowledge -based 
engineering library 12 . 

The tools 10 may also include various 
design tools, shown generally at 16, which can be 
used for this design method 20, in a manner to be 

15 described. These design tools 16 may include solid 
modeling and parametric design techniques. Solid 
modeling, for example, takes electronically stored 
vehicle design data from the vehicle platform library 
14 and standard threaded fastener data from the 

20 knowledge -based engineering library 12 and builds 
complex geometry for part-to-part or full assembly 
stress analysis. Several modeling programs are 

commercially available and generally known to those 
skilled in the art. 

25 The parametric design technique is used in 

the electronic construction within a computer system 
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22 to be described, for a geometrically defined 
vehicle system, such as the threaded fastener. As a 
particular dimension or parameter is modified for a 
particular feature of the threaded fastener or 
5 joining structure, the computer system 22 is 
instructed to regenerate a new geometric model. The 
knowledge -based engineering library 12 is used to 
control and limit the design process in accordance 
with predetermined design parameters, to be 
10 described. 

The tools 10 also include various computer- 
aided engineering (CAE) analysis tools 18. One 
example of a CAE analysis tool 18 is a finite element 
analysis software program, such as NASTRAN, or 

15 ABAQUS . Another example of a CAE analysis tool 18 is 
a visualization software program. Finite element 
analysis is used to analyze stress in a structure and 
to predict fatigue failures in areas of the structure 
with high stresses. The finite element analysis 

20 program receives as an input a data file containing a 
mesh model of the structure to be analyzed, as in 
known in the art. The output from the finite element 
analysis may be presented visually using the 
visualization software program. 

25 The tools 10 further include the computer 

system 22 as is known in the art to implement a 
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method of generating a finite element mesh for a 
threaded fastener and joining structure according to 
the present invention for use by the finite element 
analysis software program. The computer system 22 
5 includes a processor and a memory 24a, which can 
provide a display and animation of a system, such as 
the threaded fastener and joining structure, on a 
display device such as a video terminal 24b. 
Parametric selection and control for the design can 

10 be achieved by a user 26, via a user interactive 
device 24c, such as a keyboard or a mouse. The user 
2 6 inputs a set of parameters and set of instructions 
into the computer system 22 when prompted to do so by 
the method 20. The set of parameters and the set of 

15 instructions may be product specific, wherein other 
data and instructions non-specific to the product may 
already be stored in the computer system 22. 

One example of an input method is a pop-up 
window with all current parameters, including an on- 

20 line description for the parameter and a current 
value therefore. For example, parametric values may 
be chosen from a table within a two-dimensional mode, 
since some vehicle designers prefer to view an 
assembly in sections which can be laid out on a 

25 drawing . 



12 

Once the computer system 22 receives the 
set of parameters and instructions from a user 26, 
the computer system 22 utilizes a system and method, 
discussed in detail subsequently, to generate the 
5 finite element mesh. 

Advantageously, the computer implemented 
system and method of generating a finite element mesh 
for a threaded fastener and joining structure, 
according to the present invention, combines all of 

10 the foregoing to provide an efficient, flexible, and 
rapid design. Further, a data file 28 containing a 
finite element mesh of the threaded fastener and 
joining structure is an output of the system and 
method 20, and the data file 28 is available for 

15 further analysis and study. 

Referring to FIGS. 2, a threaded fastener 
and joining structure assembly 30, according to the 
present invention, is illustrated. It should be 
appreciated that the threaded fastener and joining 

20 structure assembly 3 0 are part of a vehicle (not 
shown), and in particular an automotive vehicle. It 
should also be appreciated that the threaded fastener 
32 of this example is a shoulder bolt, and the 
joining structure 34 is an aluminum end cap bolted to 

25 a steel engine block, as are both known in the art. 

The threaded fastener 32 includes a head 3 6 
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extending radially and a shaft 38 extending axially 
from the head 36. A portion of the head 3 6 near the 
shaft 38 forms a radially extending flange 40. A 
proximate end of the shaft 38 includes a shoulder 
5 portion 42, which extends radially a predetermined 
distance and axially a predetermined distance. The 
shaft 38 includes a nonthreaded portion 44 extending 
axially a predetermined distance and a threaded 
portion 46 extending axially a predetermined 

10 distance. A thread 48 on the shaft 3 8 has a profile 
with a root 50 and a peak 52 . A distal end of the 
shaft 3 8 includes a tip portion 54, which extends 
radially a predetermined distance and axially a 
predetermined distance. 

15 The joining structure 34 in this example 

includes two blocks that are secured together by the 
threaded fastener 32. Preferably a clearance hole 
block 56 is the end cap and a threaded block 58 is 
the engine block. The clearance hole block 56 

20 includes a clearance block aperture 60 extending 
therethrough the clearance hole block 56 for 
receiving the shoulder 42 and unthreaded portions 44 
of the shaft 38 of the threaded fastener 32, known in 
the art as a clearance hole. 

25 The threaded block 58 includes a threaded 

block aperture 62 extending partially therethrough 
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the threaded block 58 for receiving the threaded 
portion 46 and tip 54 of the shaft 38 of the threaded 
fastener 32. The threaded block 58 also includes a 
chamfer 64 surrounding the threaded block aperture 
5 62, to provide a counter bore. A proximate wall 
portion of the threaded block aperture 66 is 
threaded, to cooperatively engage the threaded shaft 
4 6 of the threaded fastener 32. A thread 68 in the 
^ threaded wall 66 also includes a root 7 0 and a peak 

10 72. Preferably, the threads 46, 68 have a different 
thread profile. A distal portion of the threaded 

— 

block aperture 74 is not threaded, to receive the tip 
54 of the threaded fastener 32, and includes a 

CO 

IZ cylindrically shaped portion 76 adjoining a conically 

;z 15 shaped portion 78. An end of the threaded block 58 

known as a threaded block cap 80 is solid. 

Referring to .FIG. 3 A, a flowchart of a 
system and method of generating a finite element mesh 
for a threaded fastener and a joining structure 
20 assembly 30, according to the present invention, is 
illustrated. Advantageously, the system and method 
embeds knowledge, guidelines and lessons learned from 
design, engineering and manufacturing experts to 
interactively mesh a threaded fastener 32 and joining 
25 structure 34. The methodology begins in bubble 200, 
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when it is called for by the user 26. The 
methodology advances to block 205. 

In block 205, the user 26 determines a type 
of fastener and joining structure assembly 30, for 
5 the threaded fastener 32 , clearance hole block 56 and 
threaded block 58 of this example. Advantageously, a 
portion of the joining structure 34 is used to build 
a sub-model of the overall or global joining 
structure, which in this example is the end cap and 

10 engine block, to improve the accuracy of the mesh. 
The size of the global model is large in terms of 
degrees of freedom in analysis calculations, and can 
be meshed using a coarse mesh and by simplifying the 
behavior of the threaded fastener. The sub model is 

15 more accurate since it is a more refined mesh, and 
uses displacement and thermal boundary conditions on 
the sub model from the global model analysis. The 
methodology advances to block 210. 

In block 210, the methodology generates a 

20 mesh of the model by advancing to circle A of FIG. 
3B, to be described. Preferably, a data file 

containing a mesh of the threaded fastener 32, 
clearance hole block 56 and threaded block 58 in a 
format usable for further analysis, is returned from 

25 circle A. The methodology advance to block 215. 
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In block 215, the methodology prepares an 
input file using the data file containing the mesh of 
the threaded fastener 32, clearance hole block 56 and 
threaded block 58 and additional input data from 
5 block 220. For example, the input data may include 
information regarding a material property, or a load 
or a boundary condition for the threaded fastener 32. 
The methodology advances to block 225. 

In block 225, the methodology uses the 

10 input file to evaluate stresses in the model of the 
threaded fastener 32, clearance hole block 56 and 
threaded block 58 using a computer aided engineering 
analytical tool 18, such as finite element analysis. 
In this example, a commercially available finite 

15 element analysis software program such as ABAQUS by 
Hibbett, Karlsson & Sorenen, Inc. is utilized. The 
output of the finite element analysis software 
program includes representative data concerning the 
stress, strain, contact pressures, and other state 

20 variables of the threaded fastener and joining 
structure assembly 30. The methodology advances to 
block 230. 

In block 23 0, the methodology presents the 
output of the finite element analysis to the user 2 6 
25 in a usable format, such as by using a visualization 
software program. An example of a commercially 
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available visualization software program is EnSight 
by Computational Engineering International, or IDEAS 
by Structural Dynamics Research Corporation. The 
visualization software program generates images of 
5 the threaded fastener and joining structure assembly 
3 0 to indicate stress concentrations, such as by 
contour curves of stress or strain, contact pressures 
and other state variables . The methodology advances 
to block 235. 

10 In block 235, the methodology uses the 

output of the finite element analysis to predict a 
fatigue life of the model of the threaded fastener 
and joining structure assembly 30. For example, a 
software program such as FEMFAT by Steyr-Daimler can 

15 be used to perform fatigue calculations on the stress 
fields to predict the life of the threaded fastener 
32, clearance hole block 56 or threaded block 58 
under various conditions. The methodology advances 
to diamond 240. 

20 In diamond 240, the user 2 6 determines if 

the model of the threaded fastener and joining 
structure assembly 3 0 should be modified based on the 
output of the finite element analysis or the fatigue 
analysis. For example, if the fatigue life of the 

25 threaded fastener 3 2 is too short or the stresses too 
high, a parameter describing the threaded fastener 32 
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can be modified. Similarly, if the stress fields and 
fatigue life are acceptable, another material for the 
clearance hole block 56 or threaded block 58 can be 
evaluated for weight savings purposes . 
5 If the model should be modified, the 

methodology advances to block 245. In block 245, a 
parameter describing the model is modified and the 
methodology returns to block 205 and continues. 
Returning to diamond 240, if the user 2 6 determines 

10 not to modify the model, the methodology advances to 
bubble 250 and ends. Advantageously, the methodology- 
can be used to rapidly generate a finite element 
mesh, analyze a model, and consider design variations 
without building a physical prototype. 

15 Referring to FIG. 3B, the methodology 

generates a mesh of the model beginning in circle A, 
shown at 275, and continuing to block 280. In block 
280, the user 26 specifies parameters describing the 
threaded fastener 32, and clearance hole block 56 and 

20 threaded block 58. An example of parameters for 

describing the threaded fastener 32 include a height, 
width and geometric progression factor for a mesh in 
the bolt head 36, shoulder 42, shaft 3 8 and tip 54. 
Other examples of parameters are a thread pitch, 

25 profile type, length, diameter and number of 
revolutions for thread growth. 



An example of a parameter for describing 
the threaded block 58 includes a counter bore depth, 
diameter and geometric progression factors. Another 
example of a parameter for the describing threaded 
5 block 58 includes block shape and a profile type, 
length, diameter and number of revolutions for thread 
growth of the threaded portion 46. An example of a 
parameter for describing the clearance hole block 56 
includes block shape and aperture size and location. 

10 The methodology advances to block 285. 

In block 285, the methodology advances to 
Circle B, to be described with respect to FIG. 3C, to 
create a mesh for the threaded fastener 32. The 
methodology advances to block 290. In block 290, the 

15 methodology advances to circle C, to be described 
with respect to FIG. 3D, and creates a mesh for the 
clearance hole block 56. The methodology advances to 
block 295. In block 295, the methodology advances to 
circle D, to be described in FIG. 3E, and creates a 

20 mesh for the threaded block 58. The methodology 
advances to block 300. 

In block 3 00, the methodology checks the 
data for the mesh of the threaded fastener 32, 
clearance hole block 56 and threaded block 58. For 

25 example, the mesh may be compared to a predetermined 
rule maintained in the knowledge based engineering 
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library 12 . An example of a rule is whether data is 
consistent. Another example is if all element faces 
have an approximate regular shape. Still another is 
if solid elements have a positive volume. The 
5 methodology advances to block 305. 

In block 3 05, the methodology creates an 
output file containing a mesh model of the threaded 
fastener 32 disposed within the joining structure 34. 
The output file contains the meshes for the threaded 

10 fastener 32, clearance hole block 56 and threaded 
block 58. The methodology advances to circle 310 and 
returns to the entry point. An example of a mesh 
model of a threaded fastener and joining structure 
assembly 3 0 is illustrated in FIGS. 4A through 4C . 

15 It should be appreciated that the model is presented 
in sections for visualization purposes. 

Referring to FIG. 3C, the methodology 
defines a mesh for the threaded fastener 32 by 
meshing each geometric entity of the threaded 

20 fastener 32 separately. In this example, the 
geometric entities of the threaded fastener 32 
include the head 36, flange 40, shoulder 42, shaft 
38, thread 40 and tip 54 portions respectively. The 
methodology advances to block 330. 

25 In block 330, the methodology advances to 

circle E, to be described with respect to FIG. 3F, to 



mesh the head 36. The methodology advances to block 
33 5 and advances to circle E in FIG. 3F to mesh the 
flange 40. The methodology advances to block 340 and 
advances to circle E in FIG. 3F to mesh the shoulder 
42. The methodology advances to block 345 and 
advances to circle F shown in FIG. 3G, to be 
described, to mesh the shaft 38. 

The methodology advances to block 3 50 and 
advances to circle G shown in FIG. 3H to be 
described, to mesh a cylindrical to helical 
transition between the unthreaded portion of the 
shaft 44 and the threaded portion 46 of the shaft 38. 
The methodology advances to block 3 55 and advances to 
circle H shown in FIG. 31 to be described, to mesh 
the thread 48. The methodology advances to block 360 
and advances to circle G shown in FIG. 3H, to mesh a 
helical to cylindrical transition between the 
threaded portion 46 of the shaft 3 8 and the tip 54 of 
the shaft 38. The methodology advances to block 365 
and advances to circle E shown in FIG. 3F, to mesh 
the tip 54. An example of a mesh of the threaded 
fastener 32 is illustrated in FIG. 5A. The 
methodology advances to bubble 37 0 and returns to the 
entry point . 

Referring to FIG. 3D, the methodology 
generates a mesh for the clearance hole block 56 as 
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shown in FIGS. 4A and 4B. The methodology begins in 
circle C shown at 400 and continues to block 405. 
Preferably, the clearance hole block 56 is 
represented by a cylinder, and a cylindrical meshing 
5 technique is used to generate the mesh. In block 405, 
the methodology generates nodes for the clearance 
hole block . For example , a system of nodes are 
arranged in space in the x, y and z coordinate 
directions using a combination of linear and 

10 geometric mathematical equations, as is known in the 
art. By defining a distance between the nodes in the 
x, y and z directions of the cylindrical coordinate 
system, the nodal positions for a node are determined 
in a sequential manner. The methodology advances to 

15 block 410. 

In block 410, the methodology defines 
elements by interconnecting sets of nodes to create 
an element . For example, the nodes are 

interconnected in a varying numerical sequence and 

20 identified, such as by assigning a number, to 
indicate which nodes are associated with which 
element. The methodology advances to diamond 415. 

In diamond 415, the user 2 6 determines 
whether to represent a contact surface shown at 84 in 

25 2 -dimensions . If the user 26 determines to represent 
a contact surface 84 in 2 -dimensions , the methodology 
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advances to block 42 0. In block 420 , the user 2 6 
selects areas of the model that are in contact, and 
the methodology defines the contact area 84 in two 
dimensions. Examples of contact areas 84 include an 
5 interface between the clearance hole block 56 and 
threaded block 58. Another contact area 84 is 
between a flange 40 of the threaded fastener 32 to a 
surface of the clearance hole block 56. The 
methodology advances to bubble 435 and returns. 

10 Referring back to diamond 415 , if the user 

26 determines not to represent a contact surface 84 
2-dimensionally, the methodology advances to diamond 
425. In diamond 425, the user 26 determines whether 
to represent a contact surface in 3 -dimensions . If 

15 the user 2 6 determines to represent a contact surface 
in 3 -dimensions , the methodology advances to block 
42 0 and the user 2 6 selects areas of the model that 
are in contact, and the methodology defines the 
selected contact surface 84 in three dimensions. The 

20 methodology advances to bubble 435 m and returns. 
Returning to diamond 425, if the user 26 determines 
not to represent the contact surface 84 in three 
dimensions, the methodology advances to bubble 435 
and returns to the entry point . 

25 Referring to FIG. 3E # the methodology 

creates a finite element mesh for the threaded block 
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58. The methodology begins in circle D, shown at 
450, when called for. Advantageously, the 

methodology generates a mesh for each geometric 
entity of the threaded block 58 including the counter 
5 bore, threaded section, clearance hole, clearance 
hole tip, and threaded block top as shown in FIGS. 4A 
through 4C and 6A through 6C. The methodology 
advances to block 455. In block 455, the methodology 
generates a mesh of the threaded block counter bore 

10 64. The methodology advances to circle E, to be 
described in FIG. 3F, to create the mesh of the 
threaded block counter bore 64 . The methodology 
advances to block 460. 

In block 460, the methodology generates a 

15 mesh of the cylindrical coordinate system to helical 
coordinate system transition between the counter bore 
64 in the threaded block 58 and the threads 68 in the 
threaded block 58 . The methodology advances to circle 
G shown in FIG. 3H to be described, to mesh the 

20 cylindrical to helical transition. The methodology 
advances to block 465. 

In block 465, the methodology generates a 
mesh of the threaded block thread 68. The 
methodology advances to circle I to be described in 

25 FIG. 3 J to create the mesh of the threaded block 
thread 68. The methodology advances to block 470. In 
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block 470, the methodology generates a mesh of the 
helical to cylindrical transition between the 
threaded block threads 68 and the distal portion of 
the threaded block aperture 74. The methodology 
5 advances to circle G shown in FIG. 3H to be 
described, to mesh the helical to cylindrical 
transition. The methodology advances to block 475. 

In block 475, the methodology generates a 
mesh of the cylindrical portion of the distal portion 

10 of the threaded block aperture 74. The methodology- 
advances to circle E to be described in FIG. 3F to 
create the mesh of the cylindrical portion 76 of the 
distal portion of the threaded block aperture 74. 
The methodology advances to block 480. In block 480, 

15 the methodology generates a mesh of the conical 
portion 78 of the distal portion of the threaded 
block aperture 74. The methodology advances to circle 
E to be described in FIG. 3F to create the mesh of 
the threaded block conical portion 78 of the distal 

20 portion of the aperture 74. The methodology advances 
to block 485. 

In block 485, the methodology generates a 
mesh of the cap 80 of the threaded block 58. 
Preferably, the cap 80 of the threaded block 58 is a 

25 solid structure. The methodology advances to circle 
E to be described in FIG. 3F to create the mesh of 
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the cap 80 of the threaded block 58. The methodology 
advances to bubble 490 and returns to the entry- 
point . 

Referring to FIG. 3F, the methodology 
5 generates a mesh for portions of the threaded 
fastener 3 2 and threaded block 58 that can be 
geometrically represented by a cylinder. The 
methodology begins in circle E shown at 500 and 
continues to block 505. Preferably, a cylindrical 

10 meshing technique, as described with respect to FIG. 
3D, is used to generate the mesh. In block 405, the 
methodology generates nodes in the radial, tangential 
and vertical directions for a cylindrical mesh. For 
example, a system of nodes are arranged in space in 

15 the x, y and z coordinate directions using a 
combination of linear, non- linear, step function and 
geometric progressions mathematical equations, as is 
known in the art. The methodology advances to block 
510. 

20 In block 510, the methodology generates 

elements by interconnecting sets of nodes to create 
an element. For example, the nodes are 

interconnected in a varying numerical sequence and 
identified, such as by assigning a number, to 

25 indicate which nodes are associated with which 
element. The methodology advances to diamond 515. 
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In diamond 515 # the user 26 determines 
whether to represent a contact surface 84 as 
previously described in 2 -dimensions . If the user 26 
determines to represent a contact surface 84 in 2- 
5 dimensions, the methodology advances to block 52 0. 
In block 520, the user 26 selects areas of the model 
that are in contact, and the methodology defines the 
selected contact surface 94 in two dimensions. The 
methodology advances to bubble 53 5 and returns . 

10 Referring back to diamond 515, if the user 

2 6 determines not to represent a contact surface 2- 
dimensionally , the methodology advances to diamond 
525. In diamond 525, the user 2 6 determines whether 
to represent a contact surface 84 in 3 -dimensions . 

15 If the user 26 determines to represent a contact 
surface 84 in 3 -dimensions, the methodology advances 
to block 530 and the user 26 selects areas of the 
model that are in contact, and the methodology 
defines the selected contact surface 84 in three 

20 dimensions. The methodology advances to bubble 535 
and returns. Returning to diamond 525, if the user 
2 6 determines not to represent the contact surface 84 
in three dimensions, the methodology advances to 
bubble 535 and returns to the entry point. 

25 Referring to FIG. 3G, the methodology 

generates a mesh for the unthreaded portion of the 
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shaft 44 for the threaded fastener 32. The 
methodology begins in circle F shown at 550 and 
continues to block 555. Preferably, the unthreaded 
portion of the shaft for the threaded fastener 32 is 
5 geometrically represented by a cylinder, and a 
cylindrical meshing technique is used to generate the 
mesh as described with respect to FIG. 3D. In block 
555, the methodology generates nodes for the shaft 3 8 
in the radial, tangential and vertical directions for 

10 a cylindrical mesh. For example, a system of nodes 
are arranged in space in the x, y and z coordinate 
directions using a combination of linear, non-linear, 
step function and geometric progression mathematical 
equations, as is known in the art. The methodology 

15 advances to block 560. 

In block 560, the methodology generates 
elements by interconnecting sets of nodes to create 
an element. For example, the nodes are 

interconnected sequentially and identified, such as 

20 by assigning a number, to identify which nodes are 
associated with which element. The methodology 

advances to diamond 565. 

In diamond 565, the user 2 6 determines 
whether to define a pretension section in the shaft 

25 38 that has a preload. If the user 26 determines to 
define a pretension section, the methodology advances 



29 

to block 570. In block 570, the user 26 identifies 
the pretension portion of the shaft 38. For example, 
the user 26 may specify a load on a portion of the 
shaft 38, such as the midpoint. The methodology 
5 advances to diamond 575 and continues. If the user 
2 6 determines not to define a pretension section, the 
methodology advances to diamond 575 and continues. 

In 575, the user 26 determines whether to 
represent a contact surface 84 in 2-dimensions . If 

10 the user determines to represent a contact surface 84 
in 2-dimensions, the methodology advances to block 
580 and the user 2 6 selects a surface of the model 
that are in contact, and the methodology defines the 
selected contact surface 84 in two dimensions. The 

15 methodology advances to bubble 595 and returns. 

Referring back to diamond 575 , if the user 
2 6 determines not to represent a contact surface 84 
2-dimensionally, the methodology advances to diamond 
585. In diamond 585, the user 26 determines whether 

20 to represent a contact surface 84 in 3 -dimensions . 
If the user 26 determines to represent a contact 
surface 84 in 3-dimensions , the methodology advances 
to block 590 and the user 26 selects areas of the 
model that are in contact, and the methodology 

25 defines the contact surfaces 84 in 3-dimensions. The 
methodology advances to bubble 595 and returns. 
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Returning to diamond 585, if the user 26 determines 
not to represent the contact surface 84 in three 
dimensions, the methodology advances to bubble 595 
and returns . 

5 Referring to FIG. 3H, the methodology- 

generates a mesh for the cylindrical to helical 
transition between a non-threaded portion and a first 
thread 46a of the threaded portion 46 of the threaded 
fastener 32 or threaded block 58. An example of a 

10 cylindrical to helical transition in the threaded, 
fastener 32 is illustrated in FIG.5B. It should be. 
appreciated that the helical to cylindrical 
transition is a reverse process of the cylindrical to 
helical transition, as illustrated for the threaded 

15 fastener 32 in FIG. 5E. Advantageously, the vertical 
height between the nodes, represented by DZ, changes 
as a function of the sweep around the axis in the 
vertical axis in the DY direction. As the position 
of the nodes are generated, they move upwards such 

20 that after one revolution, the last set of nodes is 
twice as high as the first set of nodes. This yields 
a part that is flat on the bottom in the cylindrical 
coordinate system and "lock washer" shaped in the 
helical coordinate system. The methodology begins in 

25 circle G shown at 600 and continues to block 605. 
Preferably, the generated mesh provides a smooth 



31 

transition from a cylindrical coordinate system to 
helical coordinate system. In block 605, the 
methodology generates nodes for the transition part 
by stretching the hexahedral elements in the vertical 
5 direction. For example, the nodes are arranged in 
space in the x, y and z coordinate directions using a 
combination of linear, non-linear, step function and 
O geometric progression mathematical equations, as is 

IP known in the art . The methodology advances to block 

10 610. 

]2 l n block 610, the methodology generates 

* = elements by interconnecting each set of nodes in a 

stretched position to create an element. For 
^2 example, the nodes are interconnected in a varying 

15 numerical sequence and identified, such as by 
assigning a number, to indicate which nodes are 
associated with which element. When the revolution is 
complete, the hexahedral elements now fill twice the 
thread pitch. The methodology advances to block 615. 
20 In block 615, wedge elements are interleaved into the 
end of the revolution to allow one half of the 
transition layers to match with the start of the 
transition and the other half to match a helical mesh 
whose elements match the original thickness of the 
25 structure. The methodology advances to diamond 620. 
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In diamond 620, the user 26 determines 
whether to represent a contact surface 84 as 
previously described in two dimensions. If the user 
26 determines to represent a contact surface 84 in 2- 
5 dimensions, the methodology advances to block 625 and 
the user 26 selects areas of the model that are in 
contact, and the methodology defines the selected 
contact surfaces in two dimensions. The methodology 
advances to bubble 640 and returns. 

10 Referring back to diamond 620, if the user 

26 determines not to represent a contact surface 84 
2-dimensionally, the methodology advances to diamond 
63 0. In diamond 63 0, the user 2 6 determines whether 
to represent a contact surface 84 in 3 -dimensions . 

15 If the user 2 6 determines to represent a contact 
surface 84 in 3 -dimensions , the methodology advances 
to block 635 and the user 26 selects areas of the 
model that are in contact, and the methodology 
defines the selected contact surface 84 in three 

20 dimensions. The methodology advances to bubble 640 
and returns. Returning to diamond 63 0, if the user 
2 6 determines not to represent the contact surface 84 
in three dimensions, the methodology advances to 
bubble 640 and returns. 

25 Referring to FIG. 31, the methodology 

creates a mesh of the thread 46 when called for. The 
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methodology begins in circle H, shown at 650 and 
continues to diamond 655. In diamond 655 , the user 
26 determines whether to generate a mesh of the 
threads using a uniform number of elements, or a non- 
5 uniform number of elements by adding elements in 
areas of known stress concentrations. If the user 26 
determines to mesh the threads with a uniform number 
of elements, the methodology advances to block 660. 

In block 660, the methodology generates a 

10 uniform mesh of the thread 46 by first creating a 
mesh of a first thread 46a of the threaded portion 46 
of the threaded fastener 32 using progressively 
larger percentages of the elements in a radial 
direction within a helical coordinate system as shown 

15 in FIG. 5C. As the nodes grow in the radial 
direction, they must rise in the vertical direction 
to maintain the helical shape to form the actual 
thread profile. The mesh is stretched in the radial 
direction as a function of height. It should be 

20 appreciated that an axial cross-section of a uniform 
mesh will demonstrate a constant number of rows and 
columns. The mesh pattern is stretched to follow the 
thread profile. Each successive column of elements 
is stretched by a smaller percentage of a full 

25 thread . 
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The methodology advances to circle J, to be 
described with respect to FIG. 3K to generate the 
mesh of the first thread 4 6a. The methodology 
advances to block 665 and generates a mesh of the 
5 thread body 46b. The methodology advances to circle 
K, to be described with respect to FIG. 3L to 
generate the mesh of the thread body 46b. The 
methodology advances to block 670 and generates a 
mesh of the last thread 46c as shown in FIG 50. The 

10 methodology advances to circle L, to be described 
with respect to FIG. 3M, to generate the mesh of the 
last thread 46c. The methodology advances to bubble 
695 and returns. 

Returning to diamond 655, if the user 2 6 

15 determines to mesh the threads with a non-uniform 
number of elements, the methodology advances to block 
675. An example of a non-uniform mesh is illustrated 
in FIGS. 6A through 6C . Advantageously, the mesh is 
denser along the surface of the thread without 

20 changing the number of . elements of the interior 
portion of the thread that form the cylindrical to 
helical transition, by making a transition zone in 
three directions simultaneously. In block 675, the 
methodology creates a growth thread profile. It 

25 should be appreciated that a thread profile is a 
vertical cross section through a thread between 
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columns of hexahedral elements. In this example, six 
rows of pentahedral elements have been added over the 
course of thread growth. The methodology advances to 
Circle P to be described in FIG. 3Q to create a 
5 growth thread profile. The methodology advances to 
block 680. 

In block 680, the methodology creates a 
mesh for the middle bolt threads 46b using a constant 
thread profile. It should be appreciated that the 

10 last growth profile for the threaded portion 4 6 is- 
the size and shape of the full thread. This profile 
is replicated a number of times corresponding to the 
number of profiles per revolution and the length of 
the thread. Hexahedral elements are created between 

15 successive profiles. A constant growth profile is 
illustrated in FIG. 6A to show a cross section of a 
full thread. This cross-section will be repeated for 
the entire length of the full size thread. The 
methodology advances to block 685. 

20 In block 685, a mesh is created for the 

last thread 46c of the threaded fastener 32. 
Preferably, the shrink thread profile corresponds 
directly with the growth profiles for the first 
thread 46a. It should be appreciated that there are 

25 the same number of growth profiles as shrink 
profiles, only the order is reversed. The primary 
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difference is that the shrink profile transition is 
from a root diameter of the threaded fastener 32 to 
the peak diameter. The methodology advances to block 
690. 

5 In block 690, the methodology spins the 

growth profiles, constant thread profiles and shrink 
profiles into a helix, representative of their final 
position. Preferably, each profile is in the X-Z 
plane. Thus, for each successive profile, the nodes 

10 of the profile are offset another step along the Z- 
axis, and rotated another step around that axis. 

Referring to FIG. 3 J, the methodology- 
creates a mesh of the threaded block thread 68 when 
called for. The methodology begins in circle I, 

15 shown at 700 and continues to diamond 705. In 
diamond 705, the user 2 6 determines whether to 
generate a mesh of the threads 68 using a uniform 
number of elements, or using a non-uniform number of 
elements by adding elements in areas of known stress 

20 concentrations, as described with respect to the 
threaded fastener 32. If the user 26 determines to 
mesh the block threads 68 with a uniform number of 
elements, the methodology advances to block 710. 

In block 710, the methodology generates a 

25 uniform mesh of the block threads 68 by first 
creating a mesh of a first thread 68a of the threaded 
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portion of the block 3 8 using gradual stretching of 
the helical elements on a manner similar to the first 
thread 46a for the threaded fastener 32. The 
methodology advances to circle M, to be described 
5 with respect to FIG. 3N to generate the mesh of the 
first thread 68a. The methodology advances to block 
715 and generates a mesh of the thread body 68b. The 
methodology advances to circle N, to be described 
with respect to FIG. 30 to generate the mesh of the 

10 thread body s 68b. The methodology advances to block 
720 and generates a mesh of the last thread 68c. The 
methodology advances to circle O, to be described 
with respect to FIG. 3P, to generate the mesh of the 
last thread 68c. The methodology advances to bubble 

15 745 and returns to the entry point. 

Returning to diamond 705, if the user 26 
determines to mesh the threads 48 with a non-uniform 
number of elements, the methodology advances to block 
725. In block 725, the methodology creates a cut 

20 thread profile for the first thread 68a in the 
threaded block 58. It should be appreciated that a 
thread profile is a vertical cross section through a 
thread between columns of hexahedral elements . The 
methodology advances to circle Q to be described in 

25 FIG. 3R to create a cut thread profile. The 
methodology advances to block 73 0. 
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In block 73 0, the methodology creates a 
constant thread profile for the thread body 68b of 
the threaded block 58. It should be appreciated that 
the last growth profile for the threaded portion is 
5 the size and shape of the full thread. This profile 
is replicated a number of times corresponding to the 
number of profiles per revolution and the length of 
the thread. Hexahedral elements are created between 
successive profiles. The methodology advances to 
10 block 73 5. 

In block 73 5, a mesh is created for the 
last thread 68c of the threaded block 58 using a 
shrink thread profile. Preferably, the shrink thread 
profile corresponds directly with the growth profiles 

15 for the first thread 68a. It should be appreciated 
that there are the same number of growth profiles as 
shrink profiles, only the order is reversed. The 
primary difference is that the shrink profile 
transitions from a root diameter of the thread 68 to 

20 the clearance hole diameter. The methodology advances 
to block 740. 

In block 740, the methodology spins the 
growth profiles, constant thread profiles and shrink 
profiles into a helix, representative of their final 

25 position for the threaded block 58. Preferably, each 
profile is created in the X-Z plane. Thus, for each 
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successive profile, the nodes of the profile are 
offset another step along the Z-axis, and rotated 
another step around that axis. 

Referring to FIG. 3K, the methodology 
5 generates a mesh for the first thread 46b of the 
threaded fastener 32, that can be geometrically 
represented by a helix as shown in FIG. 5C . It should 
be appreciated that the first thread 46a is "grown" 
gradually from the shaft 38. A proportion of the 

10 thread profile and a proportion of the difference in 
baselines are added as the nodes are laid out around 
the vertical axis. The methodology begins in circle 
J shown at 750 and continues to block 755. 
Preferably, a helical meshing technique is used to 

15 generate the mesh. 

In block 755, the methodology generates 
nodes for the first thread 46a with gradual 
stretching on the helix. Preferably, successively 
larger percentages of the full thread profile are 

20 used to shape the nodes on the outer surface and the 
nodes closer to the axis are scaled by a percentage 
of that shape. For example, the first step of a 
thread has a profile percentage of zero resulting in 
no difference in the outer nodes from the shape of a 

25 cylinder. The nodes closer to the axis are not 
changed from those for a cylinder. In successive 



40 

steps along the first thread 46a, the profile of the 
thread is scaled by a successively larger percentage 
in the radial direction to compute a position for the 
outer nodes. Nodes closer to the axis use the shape 
of the corresponding outer nodes and scale that shape 
in the radial direction by the percentage of distance 
from the axis. Preferably, the first thread 46a grows 
gradually from a flat surface to a triangular height. 
The methodology advances to block 760. 

In block 760, the methodology generates 
elements by interconnecting sets of nodes to create 
an element. For example, the nodes are 

interconnected in a varying numerical sequence and 
identified, such as by assigning a number, to 
indicate which sets of nodes are associated with 
which element . The methodology advances to diamond 
765. 

In diamond 7 65, the user 2 6 determines 
whether to represent a contact surface 84 in 2- 
dimensions. If the user 26 determines to represent a 
contact surface 84 in 2 -dimensions, the methodology 
advances to block 77 0 and the user 26 selects areas 
of the model that are in contact, and the methodology 
defines the selected contact surfaces 84 in two 
dimensions. An example of a contact surface 84 is 
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the fastener thread 46 and block thread 68. The 
methodology advances to bubble 785 and returns. 

Referring back to diamond 765 # if the user 
2 6 determines not to represent a contact surface 2- 
5 dimensionally, the methodology advances to diamond 
775. In diamond 775, the user 26 determines whether 
to represent a contact surface 84 in 3-dimensions . 
If the user 2 6 determines to represent a contact 
surface 84 in 3-dimensions, the methodology advances 

10 to block 780 and the user 26 selects areas of the 
model that are in contact, and the methodology 
defines the selected contact surfaces 84 in three 
dimensions. The methodology advances to bubble 785 
and returns. Returning to diamond 775, if the user 

15 determines not to represent the contact surface in 
three dimensions, the methodology advances to bubble 
785 and returns to the entry point. 

Referring to FIG. 3L, the methodology 
generates a mesh for the thread body 46b of the 

20 threaded fastener 32, that can be geometrically 
represented by a helix. Preferably, a helical meshing 
technique is used to generate the mesh as described 
in FIG. 3K, with constant progressive growth or 
stretching of the nodes. The methodology begins in 

25 circle K shown at 800 and continues to block 805. In 
block 805, the methodology generates nodes for the 
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thread body 46b by constant progressive growth on the 
helix by moving nodes upwards to use the helical 
coordinate system in a constant manner, as is known 
in the art. The methodology advances to block 810. 

In block 810, the methodology generates 
elements by interconnecting each of the nodes to 
create an element. For example, the nodes are 
interconnected in a varying numerical sequence and 
identified, such as by assigning a number, to 
indicate which nodes are associated with which 
element. The methodology advances to diamond 815. 

In diamond 815, the user 26 determines 
whether to represent a contact surface 84 in 2- 
dimensions. If the user 26 determines to represent a 
contact surface 84 in 2 -dimensions , the methodology 
advances to block 82 0 and the user 2 6 selects areas 
of the model that are in contact, and the methodology 
defines the selected contact surfaces 84 in two 
dimensions. The method advances to bubble 83 5 and 
returns . 

Referring back to diamond 815, if the user 
2 6 determines not to represent a contact surface 2- 
dimensionally, the methodology advances to diamond 
82 5. In diamond 825, the user 2 6 determines whether 
to represent a contact surface in 3-dimensions . if 
the user 26 determines to represent a contact surface 
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84 in 3 -dimensions, the methodology advances to block 
83 0 and the user 2 6 selects areas of the model that 
are in contact, and the methodology defines the 
selected contact surfaces 84 in three dimensions. 
5 The methodology advances to bubble 835 and returns. 
Returning to diamond 825, if the user 26 determines 
not to represent the contact surface 84 in three 
dimensions, the methodology advances to bubble 83 5 
and returns. 

10 Referring to FIG. 3M, the methodology 

generates a mesh for the last thread 46c of the 
threaded fastener 32, that can be geometrically 
represented by a helix. Preferably, a helical meshing 
technique is used to generate the mesh as described 

15 with respect to FIG. 3K, except that the last thread 
"shrinks" into the shaft 38. The methodology begins 
in circle L shown at 850 and continues to block 855. 
In block 855, the methodology generates nodes for the 
last thread 4 6c. For example, nodes are moved 

20 radially inwards toward an axis of the threaded 
fastener 32 by constant shrinking on the helical 
coordinate system, as is known in the art. 
Preferably, the last thread 46b shrinks constantly 
from a triangular height to a flat surface. The 

25 methodology advances to block 860. 
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In block 860 , the methodology generates 
elements by interconnecting sets of nodes to create 
an element . For example, the nodes are 

interconnected in a varying numerical sequence and 
5 identified, such as by assigning a number, to 
indicate which sets of nodes are associated with 
which element. The methodology advances to diamond 
865. 

In diamond 865, the user 2 6 determines 

10 whether to represent a contact surface 84 in 2- 
dimensions. If the user 26 determines to represent a 
contact surface 84 in 2 -dimensions , the methodology 
advances to block 87 0 and the user 26 selects areas 
of the model that are in contact, and the methodology 

15 defines the selected contact surfaces 84 in two 
dimensions. The methodology advances to bubble 890 
and returns . 

Referring back to diamond 865, if the user 
2 6 determines not to represent a contact surface 84 

20 2-dimensionally, the methodology advances to diamond 
875. In diamond 875, the user 26 determines whether 
to represent a contact surface in 3-dimensions . If 
the user 2 6 determines to represent a contact surface 
84 in 3-dimensions, the methodology advances to block 

25 880 and the user 26 selects areas of the model that 
are in contact, and the methodology defines the 
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contact surfaces 84 in three dimensions. The 
methodology advances to bubble 895 and returns. 
Returning to diamond 875, if the user 26 determines 
not to represent the contact surface 84 in three 
5 dimensions, the methodology advances to bubble 890 
and returns . 

Referring to FIG. 3N, the methodology 
generates a mesh for the first thread 68a of the 
threaded portion of the threaded block 58, that can 

10 be geometrically represented by a helix. Preferably, 
a helical meshing technique is used to generate the 
mesh as described with respect to FIG 3K. However, 
for the threaded block the first thread 68a resembles 
a full helical thread cut with a horizontal plane. 

15 The cut is formed by changing the thread profile 
used, such that a lower portion of the profile is cut 
away. As the first thread is swept around one 
revolution, the nodes are laid out beginning with a 
flat profile and ending with a full profile. Half 

20 way through the sweep, the lower half of the profile 
is cut away. The methodology begins in circle M shown 
at 900 and continues to block 905. In block 905, the 
methodology generates nodes for the first thread 68a 
by cutting the first thread 68a on the helix in the 

25 manner described. Preferably, the first thread 68a 
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on the block is cut into on a horizontal plane. The 
methodology advances to block 910. 

In block 910, the methodology generates 
elements by interconnecting sets of nodes to create 
5 an element. For example, the nodes are 

interconnected in a varying numerical sequence and 
identified, such as by assigning a number, to 
indicate which sets of nodes are associated with 
which element . The methodology advances to diamond 
10 915. 

In diamond 915, the user 2 6 determines 
whether to represent a contact surface 84 in 2- 
dimensions. If the user 26 determines to represent a 
contact surface 84 in 2 -dimensions , the methodology 
15 advances to block 92 0 and the user 26 selects areas 
of the model that are in contact, and the methodology 
defines the selected contact surfaces 84 in two 
dimensions. The methodology advances to bubble 935 
and returns . 

20 Referring back to diamond 915, if the user 

26 determines not to represent a contact surface 84 
2 -dimensional ly, the methodology advances to diamond 
925. In diamond 925, the user 26 determines whether 
to represent a contact surface 84 in 3 -dimensions . 

25 If the user 26 determines to represent a contact 
surface 84 in 3 -dimensions, the methodology advances 
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to block 93 0 and the user 2 6 selects areas of the 
model that are in contact, and the methodology 
defines the selected contact surfaces 84 in three 
dimensions. The methodology advances to bubble 93 5 
5 and returns . Returning to diamond 925 , if the user 
2 6 determines not to represent the contact surface 84 
in three dimensions, the methodology advances to 
bubble 935 and returns. 

Referring to FIG. 30, the methodology 

10 generates a mesh for the thread body 68b of the 
threaded block 58, that can be geometrically 
represented by a helix. Preferably, a helical meshing 
technique is used to generate the mesh as described 
with respect to the threaded fastener 32 in FIG. 3L. 

15 The methodology begins in circle N shown at 950 and 
continues to block 955. In block 955, the methodology 
generates nodes for the thread body 68b by constant 
progressive growth on the helix. For example, nodes 
are moved upwards to use the helical coordinate 

20 system in a constant progressive manner, as is known 
in the art. The methodology advances to block 960. 

In block 960, the methodology generates 
elements by interconnecting sets of nodes to create 
an element. For example, the nodes are 

25 interconnected in a varying numerical sequence and 
identified, such as by assigning a number, to 
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indicate which sets of nodes are associated with 
which element . The methodology advances to diamond 
965. 

In diamond 965, the user 2 6 determines 
5 whether to represent a contact surface 84 in 2- 
dimensions. If the user 26 determines to represent a 
contact surface 84 in 2 -dimension, the methodology- 
advances to block 970. In block 970, the user 26 
selects areas of the model that are in contact, and 

10 the methodology defines the selected contact surfaces 
84 in two dimensions. The methodology advances to 
bubble 985 and returns. 

Referring back to diamond 965, if the user 
26 determines not to represent a contact surface 84 

15 2 -dimensional ly, the methodology advances to diamond 
97 5. In diamond 975, the user 2 6 determines whether 
to represent a contact surface 84 in 3 -dimensions . 
If the user 2 6 determines to represent a contact 
surface 84 in 3 -dimensions , the methodology advances 

20 to block 980. In block 980, the user 26 selects 
areas of the model that are in contact, and the 
methodology defines the selected contact surfaces 84 
in three dimensions. The methodology advances to 
bubble 985 and returns. Returning to diamond 975, if 

25 the user 2 6 determines not to represent the contact 
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surface 84 in three dimensions, the methodology 
advances, to bubble 985 and returns. 

Referring to FIG. 3P, the methodology- 
generates a mesh for the last thread 68c of the 
5 threaded block 58, that can be geometrically- 
represented by a helix. Preferably, a helical meshing 
technique is used to generate the mesh as described 
with respect to FIG.3M. The methodology begins in 
circle O shown at 1000 and continues to block 1005. 

10 In block 1005, the methodology generates nodes for 
the last thread 68c by gradual w shr inking " on the 
helix. For example, by gradually moving nodes in a 
radial direction away from the axis of the threaded 
fastener 22 on the helix to use the coordinate system 

15 in a constant manner, as is known in the art. 
Preferably, the last thread shrinks constantly from a 
triangular height to a flat surface. The methodology 
advances to block 1010. 

In block 1010, the methodology generates 

20 elements by interconnecting sets of nodes to create 
an element. For example, the nodes are 

interconnected in a varying numerical sequence and 
identified, such as by assigning a number, to 
indicate which sets of nodes are associated with an 

25 element. The methodology advances to diamond 1015. 
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In diamond 1015, the user 2 6 determines 
whether to represent a contact surface 84 in 2- 
dimensions. If the user 26 determines to represent a 
contact surface 84 in 2 -dimensions, the methodology 
5 advances to block 1020 and the user 26 selects areas 
of the model that are in contact, and the methodology- 
defines the selected contact areas in two dimensions. 
The methodology advances to bubble 1035 and returns. 

Referring back to diamond 1015, if the user 

10 2 6 determines not to represent a contact surface 84 
2 -dimensional ly, the methodology advances to diamond 
1025. In diamond 1025, the user 26 determines whether 
to represent a contact surface in 3-dimensions . If 
the user 26 determines to represent a contact surface 

15 in 3-dimensions, the methodology advances to block 
103 0 and the user 2 6 selects areas of the model that 
are in contact, and the methodology defines the 
selected contact surfaces 84 in three dimensions. 
The methodology advances to bubble 1035 and returns. 

20 Returning to diamond 1035, if the user 26 determines 
not to represent the contact surface 84 in three 
dimensions, the methodology advances to bubble 103 5 
and returns. 

Referring to FIG. 3Q, a growth thread 

25 profile is created. A thread profile is a vertical 
cross section through a thread between columns of 
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hexahedral elements. The methodology begins in circle 
P shown at 1050 when called for and advances to block 
1055. In block 1055, the methodology meshes a 
starting flat profile. For example, a thread mesh is 
5 constructed by creating the nodes on a thread profile 
and constructing elements between matching nodes of 
two adjacent profiles. Preferably, the outer two 
layers of nodes are positioned so that the layers are 
thin to form boundary layer elements to better 

10 represent the high changes in stress near a surface 
of the threads, and the remaining nodes are smoothed 
into an optimal position. To connect smoothly to the 
shaft helix, the first thread profile and the bottom 
of the thread profiles are hexahedral elements . 

15 Advantageously, the smooth connection decreases the 
opportunity for stress analysis errors between the 
shaft 3 8 and threads 48 than would occur if a multi- 
point constraint is used. The first profile has a 
flat side which is a continuation of the cylinder. 

20 The outer edge of each additional profile bows 
outward in another percentage step to the size of the 
full profile.. Each growth profile is tilted to 
create a smooth transition from the shaft diameter to 
the thread root diameter. The first few profiles are 

25 similar to the first one, except that the outer edge 
is bowed more outward into the shape of the thread. 
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The methodology advances to block 1060 and 
meshes a transition profile. Preferably, the 

transition profile allows changing from the flat 
elements of the shaft to the more dense elements of 
the full thread. For example, the full hexahedron is 
shown at 1150 in FIG. 7 and is represented by points 
ABCDEFGH, which can be divided into smaller 
hexahedrons Aijkmnpr, iBCjnFGp, kjCDrpGH and mnprEFGH 
shown at 1155. Faces BCGF , DCGH and EFGH shown at 
1160 match with full size hexahedra that are 
adjacent. The patterns on faces ABCD, ADHE, and ABFE 
must now propagate through the mesh as smaller 
elements. Thus, large elements have transitioned to 
smaller elements in all three directions. The nodes 
mnpr will be in a position halfway between two 
profiles that would contain nodes ABCDijk in one 
profile and EFGH in another. The full transition 
includes two of these triple transition blocks, the 
second a mirror image of the first. 

The methodology advances to block 10 65 and 
creates a mesh of a wedge profile. It should be 
appreciated that the wedge elements are used to add 
additional rows of hexahedron for thread growth. For 
example, for every few profiles, another pair of 
nodes is added to the outer shape of the thread. 
This allows the mesh to maintain a proper density for 
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analytical purposes. Preferably, the number of node 
pairs is determined from the specifications and the 
pairs are added at regular intervals through the 
thread growth process . Two more nodes on the outer 
5 shape becomes another row of nodes across the thread. 
This profile is formed after the first pair of nodes 
is inserted along the outer shape along with the new 
row of elements. The quadrilateral faces of these 
elements form the row just inside the points with 

10 five edges. Advantageously, the three edge and five 
edge points are spread farther apart allowing the 
thread room to grow in an outward direction. The. 
methodology advances to 1070. 

In block 1070, the methodology meshes the 

15 hexahedron and pentahedron elements. For example, 
hexahedral elements are built in locations where two 
consecutive profiles correspond. This correspondence 
will cover the whole profile except at transition and 
wedge profiles. The transition profile requires the 

20 creation of the transition hexahedral element. Wedge 
profiles have one row of pentahedral elements at the 
new row of nodes across the profile. The methodology 
advances to block 1075 and returns. 

Referring to Fig. 3R a cut thread profile 

25 is created for the first thread 68a of the threaded 
block 58. In this example, the first thread 68a of 
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the threaded block 58 has the appearance of being 
sheared by the plane that marks the joint with the 
counter bore 64, and the first thread 68a is built 
above this plane . Preferably, the first thread 68a 
5 is built by constructing the cross-sectional profiles 
of the mesh. These profiles are the pattern of nodes 
on a vertical cross-sectional plane through the mesh 
at regular intervals around the axis of the threaded 
block 58. Solid elements will be built between these 

10 profiles using the nodes in the profiles. 

The methodology begins in circle Q shown at 
1100 when called for and advances to block 1105. In 
block 1105, the methodology creates a mesh of a 
starting threaded block profile consisting of a 

15 simple mesh of nodes in regular rows and columns. 
The first profile consists of one row of nodes and 
the second profile will have two rows as shown in 
FIG. 8A at 1180. Additional rows are added in 
succeeding profiles until the transition profile is 

20 needed. The rows of nodes are positioned to fit into 
the space between the shear plane and the full thread 
that will be built above it. Though these profiles 
will be spaced regularly around the block, they are 
created in the x-z (radial-axial ) plane . The 

25 methodology advances to block 1110. 



55 

In block 1110, the methodology generates a 
mesh of a transition profile. The transition profile 
is similar to the transition profile described with 
respect to the threaded fastener 32, but is made up 
5 of only one of the triple transition blocks instead 
of two. In addition, more rows of nodes are added to 
the profile to bring the current number of rows to 
half of the eventual total. The methodology advances 
to block 1115. 

10 In block 1115, the methodology generates a 

mesh of the first half of the thread profile to fill 
in half of the sheared thread. Each succeeding 
profile extends farther from the thread root as more 
of the full thread appears above the shear plane. 

15 Several of the profiles contain an additional column 
of nodes as compared to the previous profile, to 
properly divide the distance between the root and the 
current extension of the thread. The space between 
rows of nodes is expanded to fill in the available 

20 space between the shear plane and the full thread 
that will be built above it. The methodology 

advances to block 1120. 

In block 1120, the methodology creates a 
mesh of a second half of the thread profile affected 

25 by the shear plane. Each of these profiles starts 
with the same pattern of nodes created in the last of 
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the first half profiles. Additional rows of nodes 
are added to fill out the profile mesh configuration 
while keeping the mesh confined against the shear 
plane. As each row is added, the transition pattern 
5 of the thread crown and of the lower root is made 
more complete. Both transition patterns are 

distorted due to the shear plane compressing the mesh 
pattern. Referring to FIG. 8B, a full size thread 
profile is illustrated at 1185. The .methodology 

10 advances to block 1125. 

In block 1125, the methodology creates a 
mesh of solid elements positioned between the 
profiles . Hexahedral elements are built where two 
consecutive profiles correspond. Where there are a 

15 different number of rows in consecutive profiles, a 
row of pentahedron elements is created. Where there 
are a different number of columns in consecutive 
profiles in the first half of the profiles, a column 
of pentahedron elements is created. If a row of 

20 pentahedron elements intersects a column of 
pentahedron elements, two tetrahedron elements are 
created to make the transition. Referring to FIG. 8C, 
a cross-section of a mesh of the threaded block is 
illustrated at 1190. A profile constrained against a 

25 shear plane is illustrated at 1192, and a full size 
profile at 1194. A shrink profile is illustrated at 
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1196. The methodology advances to bubble 1130 and 
returns . 

The present invention has been described in 
an illustrative manner. It is to be understood that 
5 the terminology, which has been used, is intended to 
be in the nature of words of description rather than 
of limitation. 

Many modifications and variations of the 
present invention are possible in light of the above 
10 teachings. Therefore, within the scope of the 

appended claims, the present invention may be 
practiced other than as specifically described. 



